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The recent Quantum Hall experiments in graphene have confirmed the theoretically well-understood 
picture of the quantum Hall (QH) conductance in fermion systems with continuum Dirac spectrum. 
In this paper we take into account the lattice, and perform an exact diagonalization of the Landau 
problem on the hexagonal lattice. At very large magnetic fields the Dirac argument fails completely 
and the Hall conductance, given by the number of edge states present in the gaps of the spectrum, 
is dominated by lattice effects. As the field is lowered, the experimentally observed situation is 
recovered through a phenomenon which we call band collapse. As a corollary, for low magnetic field, 
graphene will exhibit two qualitatively different QHE's: at low filling, the QHE will be dominated 
by the "relativistic" Dirac spectrum and the Hall conductance will be odd-integer; above a certain 
filling, the QHE will be dominated by a non-relativistic spectrum, and the Hall conductance will 
span all integers, even and odd. 



The quantum Hall effect (QHE) is one of the richest 
phenomena studied in condensed matter physics. This 
effect is characterized by certain conductance proper- 
ties in two-dimensional samples i.e. the vanishing of the 
longitudinal conductance cf^x ^ along with the onset 
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of a quantized transverse conductance cjxy = ^^7^. Re- 
cently several experimental groups have produced two- 
dimensional plane films of graphite, commonly known as 
graphene sheets^ which exhibit interesting QHE behav- 
ior. 

Graphene has a theoretical history beginning with the 
study of the band structure of this planar system in^. 
From these humble beginnings it has gone on to be stud- 
ied intensely because of its Dirac structure. The bands 
can be effectively characterized by massless (2^1) d Dirac 
fermions"^. This continuum model of graphene has been 
subsequently used to study the (2 + parity anomaly^ 
and as a model system for the relativistic quantum Hall 
effect (RQHE)^'^'^. A quantum spin Hall effect has also 
been predicted in graphene^' but the intrinsic spin or- 
bit gap is probably too small to support a measurable 
phasaiiiia. 

The latter studies were based on the recent experimen- 
tal work done on the QHE in graphene by two indepen- 
dent group&^. These two groups confirm an interest- 
ing behavior in graphene in which the transverse con- 
ductance is quantized as an integer plus a half-integer 
(^xy = + ^)4e^//i, where band and spin degeneracies 
have been taken into account. Although unrelated to the 
parity anomaly, this behavior of the Hall conductance 
was in fact obvious in the seminal work of Jackiw and 
Rebbii^. On the basis of the argument for the RQHBSi2i& 
the experimental groups conclude that this is an inter- 
esting new phenomena completely explained by the rela- 
tivistic Dirac spectrum of graphene. We want to improve 
on this argument for several reasons. For very large B 



the lattice is expected to dominate the behavior of the 
Hall conductance. In this regime the Dirac argument 
cannot be valid, since, by virtue of being a continuum ar- 
gument, it ignores lattice effects and the torus structure 
of the Brillouin zone. We will see this is indeed the case, 
and the large-5 limit does not match the Dirac argument 
prediction. On the other hand, in the experimental situ- 
ation the magnetic field is weak (with respect to the unit 
quantum flux per plaquette) and the Dirac argument ap- 
plies, it is nonetheless desirable to have a description of 
the quantum Hall effect valid for both strong and weak 
magnetic fields. At low filling, we show how graphene 
evolves from a high-5 regime with non-Dirac behavior 
to a low-5 regime with Dirac behavior through a phe- 
nomenon we dub "band collapse." Two adjacent bands 
close the gap between them across the whole Brillouin 
zone and form a new band with twice the degeneracy of 
each of the initial bands. The edge structure reflects this 
degeneracy. 

We begin with a restatement of the RQHE argument 
based on the relativistic (2+l)(i Dirac spectrum. We then 
present the exact solution of the Landau problem on the 
graphene lattice. The agreement we find between numer- 
ical diagonalization and analytic calculations done with 
Hatsugai's^^ theoretical framework lead us to our conclu- 
sions and illustrate the competition between the relativis- 
tic and non-relativistic character of the band structure of 
graphene in a magnetic field. 



THE RELATIVISTIC QUANTUM HALL 
EFFECT IN GRAPHENE 



We start with the tight-binding nearest neigh- 
bor Hamiltonian for the hexagonal lattice given by 
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Semenoflii: 

H = -tJ^AA cHA)c{A + bi) + ct(A + bi)c{A) 
+/3Elct(l)c(l) - c^{A + bi)c{^+bi) (1) 

where c{A),c{A + bi) are the annihilation operators for 
sites on sublattice A and and (3 is an energy difference 
for electrons localized on the A and B sublattices. We 
will call this term the Semenoff mass. Graphene is effec- 
tively massless which is approximated by taking /3 ^ 0. 
In this limit the band structure is gapless at two inequiv- 
alent points K = -^^{^^ 2^3 ^'^^ ~ ^^^^^ ^ 
nearest-neighbor lattice constant. Around these points, 
the Hamiltonian is described by (in the ideal case mass- 
less) Dirac fermions with4i&: 

Hk = cTxkx -\- cTyky] Hk' = —cTxkx -\- (Tyky (2) 

which act on a two-spinor wavefunction describing the 
sublattices A and 5, see Fig^. There is also an overall 2- 
fold spin degeneracy which we neglect for the remainder 
of the paper. Note that parity switches A ^ B and 
K ^ K' while time reversal switches K ^ K' . The 
Semenoff term opens a gap of value m = 2/?/V^ta at 
K and — m at i^' = —K^o time reversal symmetry is 
preserved. 



where ipn{x) are harmonic oscillator eigenstates and 
are the eigenstates of Hk with energies . Notice that 
all the energy levels are paired except the n = level. 
There is a common misconception that unpaired "zero- 
modes" occur only for a massless fermion but observe 
that for m > we have ^^^q = while for m < we 
have u'^Q = 0, so such levels are unpaired even for non- 
zero mass. In the field theory formalism, the current is 
defined to be = — ^e7^^[V^Q,, V^^] and is odd w.r.t. 
charge conjugation symmetry. We find that 

{0\A0)=p=\{N.-N+)^-^ (4) 

where and N_ are the numbers of filled positive and 
negative energy Landau levels (LL). Hence the Hall con- 
ductance is 

a,y = i(iV_ - iV+) (5) 

in units of e'^/h. Due to the unpaired level, this will 
be half-integer and the position of the unpaired level de- 
pends on the sign of eB and m as in Fig[2j. 
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FIG. 2: Zero mode in the Dirac Equation. 



FIG. 1: Graphene Lattice, BZ and One dimensional lattice 
on which the Harper equation is defined. 

Now consider one Dirac fermion at the iC-point with 
mass m in magnetic field B. The Hamiltonian is H = 
c^xkx + (^y{ky — eBx) + ma^. For e5 > 0, the eigenstates 
are 



iky / ± m V^n(^ - ^0 (^)) 



(3) 



±^/an T m ^l)n-i{x - {k)) 



This analysis is correct for the fermion located around 
the i^T-point, but as mentioned before the graphene band- 
structure contains two such fermions. For the purpose of 
being well defined, we consider a small positive Semenoff 
mass m dit K which means a small negative mass at K' . 
Consider the case of eB > 0. The Hall conductance gets 
a contribution from both fermions and is zero when the 
Fermi level is in the gap — m < /i < m and odd inte- 
ger otherwise. This is then an odd integer quantum Hall 
effect as in FigP). When the gap is vanishingly small, 
m ^ the region of zero Hall conductance becomes in- 
finitely narrow. 



with 



an = \/2\eB\n + nn? 



"0 " eB 



{k±an) 



II. HARPER EQUATION FOR GRAPHENE 

We now present a different argument that reproduces 
the experimental results and is valid for both high and 
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energy spectrum, which depends only on \Aj\ is peri- 
odic with period q. We start with the one-particle states 
= i^i{k ^ (j))cl. -lO) and act on these with the 
Hamiltonian to obtain the equation = There 

are two independent amplitude equations, one for i odd 
and one for i even: 



Ajllj2j 



(8) 



where e = E/t with E the energy. There are now two 
Harper equations for the hexagonal lattice, in contrast to 
the single Harper equation for the square lattice. After 
some manipulation we find in a transfer matrix formal- 



FIG. 3: Hall conductance as a function of the chemical poten- 
tial. The red and blue are the individual conductances from 
the two Dirac cones whereas the yellow one is the total Hall 
conductance. 



low B. The solution to this problem is to carefully exam- 
ine the band structure and edge states of graphene in a 
magnetic field with rational fiux (j) = p/q. The analysis 
is based on a generalization of Hatsugai's work^^ to the 
honey-comb lattice. The energies of the bands and edge 
states are found as zeroes of certain polynomial equa- 
tions. By using general polynomial theory we are able to 
characterize the bands, find the number of band cross- 
ings, and determine the conditions for zero modes and 
edge states. By identifying the Hall conductance as the 
winding index of the edge state around the band gap, we 
find that, as the magnetic field is decreased, the wind- 
ing number of the edge states starts taking odd-integer 
values due to electron bands collapsing in pairs. The 
theoretical spectrum is obtained, and in addition, exact 
diagonalization results are presented to support it. It 
will be evident from the calculated bandstructure that 
for large magnetic fields the Dirac argument does not ap- 
ply because the Hall conductances of bands at low-filling 
do not form a sequence of odd-integers in this case, as 
predicted by the relativistic argument. 

We use the Landau gauge Ay = Bx^ = 0, with 
B = 2^/3V3a^^ where ^ = p/q is the flux per plaquette 
(hexagon) and p, q are relatively prime integers. With 
a Peierls substitution the effect of the magnetic field is 
cjcj c\cj exp p. Adf. In this gauge ky = k is a. good 
quantum number and the Hamiltonian for each k is: 

H{k) = -^X14,2j-i^fc'2iA^-(A:)+c^^.c/e,2j+i+/i.c. (6) 
j 

where 



(7) 



Note that we have not included a mass term in our 
tight-binding Hamiltonian because graphene is essen- 
tially massless. Since Aj+q = {—l)^Aj the Hamilto- 
nian is periodic with period 2g, (Aj+2g = Aj) but the 



V^2j+l 



with 
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V'2i- 



(9) 



(10) 



As opposed to the transfer matrix for the square lat- 
tice, which hops by one site and is linear in energy 
the graphene transfer matrix hops by two sites and is 
quadratic in energy. This refiects the lattice periodicity. 
Since Mj-^q = Mj^ the periodicity of the energy spec- 
trum is q. We can now define the transfer matrix over 
the magnetic unit cell: 



M(e) 



_fMn{e) Mi2(e) 
M2i(e) M22(e) 



1...M1. (11) 



By induction we find that Mu is a polynomial of order 
(e^)^, M12 and M21 are of the form e x (e^)^"\ while M22 
is a polynomial of order (e^)*^"^. These polynomials have 
coefficients which depend on k and the magnetic ffux. 
We pick our sample of order Ly = 2ql^ commensurate 
with the magnetic unit cell, where / is a large integer and 
the factor of 2 is added because we will require periodic 
conditions il)Ly = V^o, hence Ly = {) = even. The transfer 
matrix across the length of this sample is MK From 
Hatsuga:^^ we know that the important polynomial to 
consider is: 



[M^]2i(e)=0 



(12) 



The entire spectrum of energy levels for each k value 
comes from the zeroes of this polynomial of which there 
are Lx — l. Some of these states are bulk states and others 
are edge states. We will now characterize the edge and 
bulk states (bands). 

It is easy to find one solution to Eq.^21- Simply take 
M2i{e) = and this will imply that Eq.^is satisfied 
since all upper-triangular matrices remain so when mul- 
tiplied by another upper-triangular matrix. Hatsugai 
argues that the energies of the edge states are given by 
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the zeroes of exactly this polynomial: M2i(e) = 0. Since 
M2i(e) --ex (a(e^)^~^ + 6(e^)^"^ + • • 0^ 'twere is always 
one e = solution (zero mode edge state) which does 
not disperse and 2{q — l) non-zero energy solutions (edge- 
states) which come in pairs: —fiq-i < —jiiq-2 ^ ••• < 
— Ml < < /ii < .... < /jjq-2 < l-iq-i' Depending on 
whether Mn^jii) /\Aq....Ai \ is <,>, or = 1 the edge state 
will be localized on the left edge, right edge, or be degen- 
erate with the bulk i.e. touching a bulk statei^. 

The bulk states are obtained from the lattice period- 
icity j ^ j -\- q and the Bloch condition: 



i^2q 



(13) 



with p(e) a pure imaginary phase, ie., \p{e)\ = 1. We 
also note that we have the transfer matrix equation 



^2q+l 
^2q 



1 



AqAq-i...Ai 



M 



(14) 



Therefore, combining these two, p(e) is an eigenvalue of 
the 2x2 transfer matrix 



[TrM±J{TrM)^-4\Aq....Ai\^]. (15) 



2Aq..A, 

where we have used Det[M] Det[Mq]...Det[Mi] 
\Aq...Ai\'^. It is easy to see that the Bloch condition 
|p(e)p = 1 is satisfied for (TrM)^ - 4\Aq...Ai\'^ < 0, 
based on the fact that p~^p~ = 1. Since Mn and M22 
are both polynomials of order q in the solutions are 
again paired. Let us rewrite 



2q 



{TrM{e')f-4\Aq...A^\' = ll{e'-X,), 



(16) 



with < Ai < A2 < • • • < X2q- The energy bands are 
thus 



A2j+1 < e < ^2j-\-2 

\2j < < A2 



A2j + 1 



bulk state 



gap region 



(17) 



for j = 0, 1, • • • , g — 1 and Aq = 0. The edge states lie 
in the gap region of the bulk band structure and the /i's 
are given by 



Pj ^ [A2j,A2j+i] j = I,-- - - 1. 



(18) 



We hence have 2q energy bands bounded by Aq A's, there 
are 2q — 1 gaps and 2q — 1 edge states as in Figg). 

Besides the above results of band structure, many more 
details about the spectrum can already be learned from 
the behavior of the function \AqAq_i...Ai\{k), 

1. The Hall conductance can be determined from the 
number of /c's that satisfy 

Mn(/i[,/2])M22(/i[,/2]) = \AqAq_i...Ai\\ki). 




FIG. 4: Schematic plot of the bulk band structure and edge 
states obtained from the transfer matrix formalism. Edge 
states are solid lines while bulk bands are denoted by the 
shaded areas bounded by dash and dotted lines. 



2. The first bulk eigenvalue touches the zero energy 
edge state at the k points where \Aq...Ai\{k) = 1. 

3. Bulk band width vanishes at k if \Aq...Ai\{k) = 0. 

For graphene, \AqAq-i...Ai\{k) can be explicitly writ- 
ten as 



\AqAq_,...Ar\{k) = 2^l[ 



k ^ (5 - 6j)p7r \ 



6q J 



.(19) 



The periodicity of this function is 27r/q. Hence the num- 
ber of /c's at which \AqAq-i...Ai\{k) = is equal to q 
while the number of /c's at which \AqAq_i...Ai\{k) = 1 is 
equal to 2q. 

We shall now show how to obtain the details of the 
band structure from \AqAq-i...Ai\{k). Let us assume 
that the edge state P[j/2] (k) touches the bulk at some 
point k = ki 



P[j/2]{h) = ej{ki) = ±JXj{ki) 



(20) 



where [j/2] represents the largest integer less than or 
equal to j/2 and j = l,...,2g'. Since ej{ki) is on the 
band edge, we have 

Mii(/i[,/2]) + M22(/i[,/2]) = ±2|A,A,_i...Ai|(/ci). (21) 

From the edge state condition M2i{jii[j/2]{ki)) = 0, we 
also know 

Mii(/i[,/2])M22(/i[,/2]) = \AqAq_i...AiWki). (22) 

We hence have 
Mn(/i[,/2](A:i)) = M22{p[j/2]{h)) = ±\Aq...A^\{ki) {23) 
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when the edge state touches the bulk state. Thus, we 
can determine how many times the edge state starts from 
^2[j/2] at ki and goes up in energy to touch X2[j/2]-\-i at 
some /c2 then comes down again to touch X2[j/2] at some 
ks etc. This defines the number of wrappings around the 
gap and represents the Hah conductance^^- 

As a function of the momentum k the first bulk eigen- 
value Ai might touch the zero energy line (the zero mode) 
when ei{k) = 0. This happens when 

Mn(0) = (24) 
M22(0) = (-1)^; (25) 
M2i(0) = 0. (26) 

But from the previous analysis we know that when a bulk 
state touches an edge state M22 = Hence the 

first bulk eigenvalue touches the zero energy edge state 
in 2q points in the first Brillouin zone, namely where 
|74g...Ai| = 1. This result is confirmed by our exact di- 
agonalization, which will be presented later. 

Using polynomial theory we can in fact prove a more 
stringent constraint. We separate the polynomial of 
order 2q: {TrM{e^)f - 4|A^...Aip = {TrM{e^) - 
2\Aq...Ai\){TrM{e^) + 2\Aq...Ai\). Now, denote the 
eigenvalues of the two subfactors as g for "green" and 
b for "blue" respectively(for purpose of making the con- 



nection with the plots) and put them in ascending order: 

TrM(62) -2|A,...Ai| = n,Li(e^-Af), 
TrM(a +2|A,...Ai| = nj=i(e^-A5), 

where Af < < ... < Ag and Aj < A^ < ... < A^. (For 
q = 2 the < sign changes into < due to the fact that in 
that case the system doesn't break T and we can have 
gapless states.) Depending on whether q is even or odd 
we have the following order: 

q odd : Aj < Af < Af < A^ < ... < A^ < Af 
q even : Af < A^ < A^ < Af < ... < A^ < Ag (27) 

We can see that bulk states are between "blue" and 
"green" eigenvalues |[A^, Af]|, whereas the gaps are in be- 
tween the consecutive "blue" -"blue" and "green" -"green" 
eigenvalues |[A^, A^^^^ll I [^i , A-'^;^] I- such, the width 
of a band is | \/A|^ — | • The band will become in- 
finitely thin when A^ = A^, or when |74g...74i| = 0. This 
happens at q points in the first Brillouin zone. 

As an example, we can see everything above explicitly 
for the case g = 3,p = 1. There are 2q — l = 5 edge states 
with energies — /i2(^), — Mi(^), 0, /ii(/c), /i2(^), where 



COs{k- y) 



•cos(/c- J)Ty^ 



cos {2k - — -) + cos {k - --)(2 + cos {k - — -)) + (cos {k - — -))2 
9 9 9 9 
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There are q = 3 points in the Brillouin zone where each 
band becomes infinitely thin given by 



I . . . I / . N . 47r IOtt IGtt , , 

\AsA2Ai\{k) = at = — , — , — . (28) 



The bands closest to zero energy touch the zero energy 
mode at 2g = 6 places in the Brillouin zone where 



\AsA2Ai\{k) = 1 at /c= -,— . (29) 



We also find that the condition Mii(/ii(/c)) = 
M22{/^i{k)) is satisfied at two points in the Brillouin zone, 
which means that, in the first gap, the edge state touches 
the lower band Af once and the upper band Af also once, 
hence the Hall conductance is one. This is the same for 
when the Fermi level rests in the second gap, the con- 
dition Mii(/ii(/c)) = M22{l^i{k)) being satisfied for two 
points in the Brillouin zone as well (see FigjS]). 



III. HALL CONDUCTANCE IN GRAPHENE 

This section contains the theoretical results from the 
transfer matrix approach, as illustrated in the previous 
section, and the numerical results from exact diagonal- 
ization. The Hall conductance in graphene is defined, 
as usual, as the number of times the edge state wraps 
around the gap between neighboring energy bands. The 
number of left or right edge states that traverse the en- 
tire way across the gap is the Hall conductance. We then 
look at the evolution of the bands and edge states as the 
magnetic field is varied from very strong to weak. We 
will see how the edge states and band configuration for 
strong magnetic field, which do not match experiment, 
evolve into the weak-field limit, which does match the 
experiments. 

This is accomplished in two ways: first, the "theoreti- 
cal" edge states and band structure are found by numer- 
ically solving for the zeroes of the characteristic polyno- 
mials M2i(/ii) = and (TrM(Af^)) ± 2\Aq...Ai\ = 
introduced in the previous section. We plot only the 
e > states, the negative energy states being a mirror 



image. We also confirm the theoretical picture by exact 
numerical diagonalization of the Hamiltonian matrix for 
a relatively large number of lattice sites. 

We start with q = 3 in FiglSJ. We see that the Hall 
conductance is unity for a Fermi level in either the first 
or second gap, clearly in contradiction with the Dirac 
argument which would give an = 1 or 3 depending on 
which gap. The number of bands is 2g = 6, there are 
2q — 1 = 5 gaps and edge states, q = 3 spots where each 
band becomes infinitely thin, and 2q = 6 points where 
the first band touches the zero energy mode. 




FIG. 5: Left: Theoretical edge state and band structure con- 
figuration. The edge states are in red while the bulk bands 
are in between the consecutive blue and green lines (A^A^). 
The number of times an edge state wraps around the bulk 
is the Hall conductance, which in this case is unity 1 for the 
both the first gap and the second gap. Right: The band struc- 
ture obtained from direct diagonalization, upper right just the 
e > lower right is the full spectrum which is just a mirror 
image of the e > spectrum. From now on, we will plot only 
the positive energy part. 

We now continue by decreasing the magnetic field to 
q = 4 and then g = 6, see Fig|S]. For g = 4 we have 
axy = 1 for a Fermi level in the first gap, axy = 2 for 
Fermi level in second gap, and axy = 1 for Fermi level 
in third gap. For q = 6 the sequence is axy = 1, 4, 3, 2, 1 
from the first to the fifth gap. This does not match the 
experimental observation of axy = 1, 3, 5, 7, etc. 

One crucial observation to notice is that, as we increase 
q (decrease the magnetic field), the second and third bulk 
bands become closer and closer together in energy; the 
gap between them becomes smaller and smaller over the 
whole Brillouin zone. Eventually the second and third 
bands move entirely together upon increasing q (lower- 
ing B). For g > 12, one cannot distinguish between the 
second and third band (nor can one distinguish the edge 
states between these bands). The second and third band 
have "collapsed" into a new band, a process which we call 
"band collapse." After these bands have collapsed there 
are distinguishable gaps between the first band and the 
combined band, and then between the combined band 




FIG. 6: Left: Theoretical edge state and band structure con- 
figuration for q = 4 and ^ = 6. Right: Direct diagonalization. 



and the fourth band. There are edge states between the 
top of the combined band and the fourth band, and these 
give axy = 3, for the Fermi level in what is now the sec- 
ond gap. If we then go to the next gap, this again does 
not match the experiment, with Hall conductance being 
8. 

By increasing q even further, we see that the fourth 
and fifth bands collapse in a similar fashion, and the gap 
between them vanishes uniformly across the k spectrum 
as they become a single new band. This happens around 
q = 22. The edge states between the collapsed second 
and third bands and the collapsed fourth and fifth bands 
remain the same as before, giving axy = 3 but now the 
edge states between the collapsed fourth and fifth bands 
and the sixth band give axy = 5. This process repeats 
itself while q is increased. The total number of bands 
increases when q is increased. But some of these bands 
collapse together so that we cannot distinguish them un- 
less we have infinite resolution. We present the results 
for q = 31 FiglHl. 

Upon increasing q the band collapse leads to double 
degeneracy of each of the bands except the zero energy 
band, and this gives the odd integer Hall conductance 
in graphene. This is beautifully seen as the number of 
positive or negative-slope edge states that disperse in the 
resolvable gaps. Hence the experimental situation is the- 
oretically confirmed as the weak-field limit of graphene. 

The theoretical band structure can actually be contin- 
ued to large q^ as the polynomials are well behaved. We 
give the q = 49 plot as well, where we can see all the 
odd-integer quantum hall effects from 1, 3, 5, 7, 9, 11, see 
Figinj. 

By examining the common properties of each band- 
structure plot it appears that the spectrum of bands and 
edge-states can be classified into two parts: a relativis- 
tic section and a non-relativistic section. This structure 
originates from the original tight-binding dispersion re- 
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FIG. 7: Left: Theoretical edge state and band structure con- 
figuration for q = 12. Right: Direct diagonalization. 



FIG. 9: Left: Theoretical edge state and band structure con- 
figuration for ^ = 49. 




6* and 7'^ band 
collapsed Into new 
band 



4t*i and 5* band 
collapsed into new 
band 



2"'^ and 3^^ band collapsed into new 
band 



1=t band, close to zero energy 




FIG. 8: Left: Theoretical edge state and band structure con- 
figuration for ^ = 31. 



lations without the B field, 



the first gap from the band ceiling and increases by one 
for each Landau level below the top. A similar thing oc- 
curs for the non-relativistic levels near the bottom of the 
set of bands. The crossover region is at £^ ~ where the 
band collapse occurs. 

The odd-integer sequence shown in FigjUJ is clearly 
represented in the experimental data which, as stated 
before, is in the low magnetic field limit of graphene. 
With a flux (j) = 1/q in each unit cell the magnetic field 
is -'-•^^-'-^ Tesla which is a very large magnetic field. 
For experimentally realizable magnetic fields we would 
expect q ~ 1000 and the odd-integer sequence would 
be continued to larger values. Abnormalities in this se- 
quence would not arise until more Landau levels were 
filled. Overall, there will be a sequence (possibly very 
long) of odd-integer quantum Hall conductances followed 
by conductances which do not follow a certain pattern. 
Then there will be a relativistic-non-relativistic crossover 
region where the Landau level spacings change character 
from n^/^ to n. The non-relativistic energy levels will 
then persist to higher energies. 




A 2 ^7 

4 cos^ — k 



4 cos — kx cos -ky 
2 2 ^ 



The Dirac nodes are located at 

.27r 



(30) 

(±^,0) and 



The linearized dispersion relations persist 
up to around E ^ t. Above this energy scale the bands 
become parabolic. Accordingly, in Fig. |2l (Jxy = 1,3,5 
... at low energy and the energy of bulk levels goes as 
En ~ ^/n, a feature of relativistic Landau levels. On the 
other hand axy = 1, 2, 3, .... starting from the top of the 
bands (where parabolic bandstructure is expected) and 
there is almost equal spacing between each of these Lan- 
dau levels, which is a feature of the harmonic-oscillator- 
like non-relativistic Landau levels. A axy of 1 is seen in 



A. Effect of Disorder 

We have considered the stability of the edge-states un- 
der disorder. Although not tractable analytically, we 
were able to use numerical diagonalization (which up 
to now has remarkably matched the analytic results) to 
study the introduction of disorder into the system. The 
disorder term we added to the system is: 



Hdis{k) = XdSij 



(31) 



where is a random variable with gaussian distribu- 
tion and mean 0, and I^J= 1, 2, ... L — 2, L — 1. We 
also tested a uniform distribution for the Xd with es- 
sentially the same results. For relatively high disorder 
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FIG. 10: Numerical calculation for edge states and band 
structure for q= 10 and disorder variance 0.15. 

e.g. the variance of ~ 0.15 the structure of the 
lowest energy edge state is robust (see Fig^Ol)- How- 
ever, the edge states representing higher plateaus, such 
as n = 3, 5, 7, 9, 11, . . . are washed out. Note that the 
hopping parameter is defined to be 1, so 0.15 is very 
high disorder. For lower disorder, with the variance of 
Xd ~ 0.01, all of the edge states are clearly visible up 
to the relativistic-non-relativistic crossover (see FigfTT]). 
just as in the disorder- free plots given above; e.g. as in 
Figs. 



FIG. 12: (color online) (a) Theoretical edge state and band 
structure configuration for q = 17, with Semenoff mass m — 1. 

The new Harper's equations are: 

(e + m)V^2j-i + ^j^2j + ^2j-2 = 

^jV^2j-i + (e - m)ilj2j + V^2j+i = (33) 

and the transfer matrix now becomes: 

M, = ('"-f-Y^ (34) 

As we can see, MqMq-i....Mi = (e + m) x p(^-^)(e^) 
where P^^-^\e^) is a polynomial of order g — 1 in 
e^. Hence the former zero energy edge state has now 
moved to /io = —tti. There are no edge states between 
[— |m|,+|m|], but the rest of the analysis applies. We 
plot the band structure for m = 1, = 17 (see FigfT^). 



IV. SPIN AND VALLEY SPLITTING IN THE 
n = LANDAU LEVEL 



FIG. 11: Numerical calculation for edge states and band 
structure for = 50 and disorder variance 0.01. 



B. Non-zero SemenofF term 

The previous formalism can be easily extended to in- 
corporate the case of a non-zero Semenoff mass m. As 
an example, for Boron Nitride (BN) the hamiltonian has 
the form: 

H = -t^{{clj_^C2jAj + 4iC2i+i + h.c.) 

3 

+^(4i-lC2i-l - cl.C2j)) (32) 



We now focus on the breaking of spin and/or valley 
degeneracy in the n = Landau Level. The idea of spin 
splitting is very natural since ^ ~ 2 in graphene and 
there is a large magnetic field applied perpendicular to 
the sample. Splitting the valleys however, is more subtle 
since there is no natural alternating sublattice potential 
or applied strain. We investigate the changes to the QHE 
plateau structure and edge states when these splittings 
can be resolved energetically. 

First we consider the case of only spin splitting. Due 
to the Zeeman effect the spin states in each Landau level 
will be split by g/J^sB. For the n = Landau level one 
spin state is pushed above zero energy and the other is 
pushed below zero energy. When the chemical potential 
lies in the gap at zero energy between the split spin states 
there is an additional QH plateau with axy = 0. The pic- 
ture is not quite this simple because this gap, unlike the 
Semenoff mass gap discussed above, contains edge states 
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+1 





















FIG. 13: Landau level bands and edge states around E=0 for 
(a) Spin splitting only (b) Valley splitting only (c) Spin and 
valley splitting. Quantum Hall conductances for particular 
gaps are noted by the integer labels and are in units of /h. 



which can be seen in Fig |13(a)l -'^^. Usually the presence of 
edge states in the gap signals a non-zero QH conductance 
but here there is actually one electron edge state and one 
hole edge state. These two edge states combine together 
to give zero Hall conductance but produce a non-zero 
spin-Hall conductivity since they are spin-polarized in 
opposite directions: 



(35) 



This spin current can be observed in a 4-terminal geom- 
etry or in a system with magnetic leads. 



The case where only the valleys are split in the n = 
level, no matter by what means, is very similar to the 
case of a non-zero Semenoff mass given above. As in that 
case there is a gap at zero energy leading to an additional 
zero conductance plateau, however here there are no edge 
states in the gap, thus no spin Hall conductivity. The 
band picture and sequence of quantum Hall conductances 
can be seen in Fig |13(b)] . 

Finally we come to the case where there are both spin 
and valley splittings. Gaps will appear when the n = 
level is unfilled, l/4-filled,l/2-filled, 3/4-filled, and com- 
pletely filled yielding a sequence of QH conductances 
o'xy = —2,-1,0,1,2 in units of e'^/h when the chemi- 
cal potential lies in each of these gaps. The band picture 
with each of these conductances can be seen in Fig |13(c)]. 
This sequence matches the data recently produced by^^* 
in very high magnetic fields. In the graphene sample there 
will be some small valley splitting due to imperfections 
(shear strain, impurities, or surface roughness) but not 
enough to produce a gap large enough to exhibit the 
quantum Hall effect. Since there is no applied strain we 
must look for many-body effects that would give rise to 
this splitting. 

The idea of exchange ferromagnetism, which applies in 
the non-relativistic quantum Hall effect is also applicable 
here with some differences. In a normal QH system we 
expect that for the lowest Landau level we should see 
valley-polarized ground states^^'^^. The dependence 
in the graphene Landau level spectrum should not be 
important as long as the Landau gap is large i.e. hcOc » 

/Ib^^/t. The excitation energy of skyrmions has been 
calculated in^^. However, if we considered higher Landau 
levels there would be some quantitative corrections. 

The second thing to consider is the correlation between 
the valley index and the sublattice index. For the n = 
level if an electron is in a particular valley then its spa- 
tial wavefunction resides on a single sublattice, A or B. 
If this Landau level is 1/4-filled or 3/4-filled there will 
be a valley and spin polarized ground state. The spin 
polarization is from the Zeeman splitting and the system 
will form a valley-polarized "ferromagnet"-like state due 
to exchange correlations. In this level the valley and sub- 
lattice are correlated, but they are correlated such that 
if the electrons reside in only one valley then they re- 
side on a single sublattice which minimizes the Coulomb 
interaction. This leads to a spin-polarized charge mod- 
ulation where there will be an excess of charge on one 
sublattice. This will form a weak charge density wave 
with charge density modulation where the percentage of 
charge modulation is proportional to Nq/Nt^ the amount 
of electrons in the n = Landau level divided by the to- 
tal number of electrons in the system. The electrons that 
participate in the charge modulation are effectively the 
difference between the number of electrons at half-filling 
and the number of electrons currently in the system. 

This valley polarized ground state will produce an in- 
teraction gap characterized by the energy to produce a 
charged excitation. Since there is no applied strain we ex- 
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pect that SU (2) valley skyrmions will be cheaper to cre- 
ate than particle- hole excitations^^. We do not expect to 
see full SU (4) skyrmions because the ^-factor in graphene 
is not small^^. This raises the possibility of measuring 
valley skyrmions in graphene as was recently done in 
AlGaAs^^. Since we are projecting into the n = Lan- 
dau level we can use the calculation of^^-^^ to estimate the 
spin stiffness and thus give an estima te o f the energy to 
create a skyrmion: Egk = ^^Ps = |\/7r/2(e^/e£B)- If we 
compare the energy width of the plateau of the spin-split 
states to that of the valley-split states shown in^^ they 
are roughly of the same order of magnitude. However 
Esk/igpsB) ~ 54 at 5 = 45T so our skyrmion energy is 
clearly an overestimate. For the valley skyrmions mea- 
sured in AlGaAs^^ the data also clearly shows that Esk 
is an overestimate by a factor of ~ 40 for their systems at 
zero applied strain. This factor compensates for the over- 
estimation and brings the skyrmion energy to the right 
order of magnitude. Another interesting fact is that at 
low magnetic field this valley splitting gap vanishes and 
the axy = ±1 plateaus disappear. This could be the 
result of there being two few electrons in the n = Lan- 
dau level to produce this well correlated effect. Overall 
the valley degeneracy splitting suggests that small spin- 
polarized charge density modulation or valley skyrmions 
could be measured in graphene. 

V. CONCLUSION 

We have shown that the "relativistic" quantum Hall ef- 
fect in graphene has its origin in a band-collapse picture 



where two bands become degenerate upon decreasing the 
fiux per plaquette. A series of exact results for the hon- 
eycomb lattice are given, as well as an index theorem for 
the number of dirac modes in a magnetic field. At large 
magnetic fields, the system has a transition between "rel- 
ativistic" and non-relativistic QHE. When the spin-gap 
is resolved, the system exhibits a spin-Hall effect due to 
existence of opposite spin electron and hole edge states in 
the gap. We discussed the effects of disorder and adding 
a Semenoff mass term. We concluded with discussion on 
spin and valley splitting in the n = Landau level and 
its implications for the quantum Hall effect. 



VI. NOTE 

During the preparation of this paper, we have no- 
ticed a series of other papers that have indepen- 
dently reached some of the conclusions presented in this 
manuscriptiSi24i2Si2Si2£. 
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